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Abstract. One of the key ingredients to understand the properties of neutrons starfl (NS) is the equation of state at finite 
densities far beyond nuclear saturation. Investigating the phase structure of quark matter that might be realized in the core of 
NS inspires theory and observation. We discuss recent results of our work to point out our view on challenges and possibilities 
in this evolving field by means of a few examples. 
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INTRODUCTION 

Not only due to their compactness with central densities multiples beyond nuclear saturation NS are considered to 
belong to the most interesting systems nature offers to study the properties of ultra-dense matter But it is this feature 
that provides a basis for the wide range of unique possibilities that NS offer to gain insight into the very nature 
of fundamental interactions. A fascinating and yet undecided question is whether the QCD phase transition from 
nucleonic to the fundamental quark matter (QM) degrees of freedom is realized in the core of NS. Answering this 
question in a strict sense would require the theorist to describe nuclear matter at finite densities in the framework of 
QCD, which is challenging already in vacuum. Solving QCD in medium at this fundamental level is a yet unsolved 
task even though promising work into this direction has been done. In a common approach to bypass this problem 
quark and nuclear matter EoS are modeled separately and are combined afterwards by constructing a phase transition 
thermodynamically. However, already pure quark or nuclear matter EoS are barely constrained due to the few reliable 
data at densities beyond saturation. 

This paper embraces aspects of our recent work on ultra-dense matter It starts with a discussion of potential 
promising constraints on the EoS of ultra-dense matter from NS mass and mass-radius measurements and the elliptic 
flow in heavy ion collisions. These, other constraints and more detailed discussions can be found in l2ll3l|4[]. The two 
following sections consider the EoS in nuclear and quark matter Here we focus on the importance of constraining 
in particular the nuclear EoS |2] and discuss present approaches on describing quark matter in medium. While 
deconfinement and dynamical chiral symmetry breaking are accessible by solving QCD's gap equation within the 
framework of the Dyson-Schwinger formalism liSJ, effective models of the NJL-type are a practical tool to explore the 
phase structure of quark matter at finite densities ||a,|711- Finally, we discuss a recently introduced approach to construct 
the phase transition from nuclear to quark matter which results in a possibly mixed phase of one-flavor quark matter 
on a background of nuclear matter 17|]. 



CONSTRAINTS ON THE EQUATION OF STATE 

Far from the basic perception of a NS as an object as heavy, and heavier, as our sun, bound by gravity into a sphere with 
a diameter of about 25 km and spinning around its axis up to 1000 times per second (which is already fascinating) 
NS are complex and inherently versatile objects. This is because NS have an internal structure and evolve in time. 
These two elements - structure and temporal history - are both a challenge and the key for a detailed understanding of 



' We use the generic term neutron stars for simplicity. In the context of this paper this does not exclude any exotic form of strongly interacting 
matter in the NS core such as color superconducting quark matter (l||. 



the physics of matter in NS. With increasing observational efforts especially over the last decade, which will increase 
further over the next, NS data become more and more valuable as constraints for the modeling of ultra-dense matter. 
Of particular interest are NS which significantly deviate from the norm since those define the natural limits one has to 
account for in any reliable theoretical description. 

Facing the present potential uncertainties in observation and modeling, it makes sense to start from an 'extremal set' 
of observables (largest and smallest known NS mass, radius, moment of inertia, rotational period, etc.) and to examine 
the consistency of this set, e.g. for similar underlying model assumptions. The benefit of this approach is apparent in 
two ways. First, a reliable EoS should be in agreement with all trusted NS data available. Second, taking into account 
extremal data is a valuable additional strategy to explore the real limits of NS properties and the EoS. Therefore, at 
present, the challenge is to identify both reliable "extreme" NS data and possible contradicting or complementing data 
that might narrow the limits ofNS observables and the EoS. 



Neutron Star Masses 

For the purpose of this work we focus on NS configurations which are extreme in terms of their mass and/or 
radius. Due to observational uncertainties one should be careful not to mistake these extreme values as the final 
truth on the subject. A relativistic analysis of timing observations of the pulsar PSR J 1903+0327 finds its mass to be 
1.74±O4M0. This is the largest pulsar mass that has ever been precisely measured Jsl]. Evidence for more massive 
NS exists but suffers from larger uncertainties. We consider a maximum NS mass of about 2.0 M^,, which is in 
agreement with timing measurements for PSR B1516+02B, located in the globular cluster M5. The derived pulsar 
mass is (2.08 ±0.19) M© (at 68% probabiUty) with a 95% probability that the mass of this object is above 1.72 Mq 
which deviates considerably from the average masses of binary radio pulsars, Mbrp — 1.35 ±0.04 M© ifioll . 
However, looking at the recent history and the prominent example of PSR J0751+1807 we are aware of the fact that 
the identification of NS masses is a sensitive task and might give flawed results. In the latter case the mass initially 
determined to be 2.1 ± O.2M0(1(7) had to be coiTected to a value of 1.26 ± O.2M0(la) due to new data analysis. 
There are other observations pointing towards high NS masses (2. 10 ± 0.28 Mq for EXO 0748-676,2.0 ± 0. 1 Mq for 
4U 1636-536) but for the given reason they should be treated with caution. As we will point out, the analysis of the 
elliptic flow in heavy ion collisions suggests that 2 Mq is an approximate upper limit on NS masses, but presently 
no other criterion is known, that would decisively rule out the existence of high massive NS. It will depend on future 
observations to affirmatively determine the actual upper limit on NS masses. 



Radii and Mass-Radius Relations 

Aside from NS masses, kilohertz quasi-periodic brightness oscillations (QPOs) seen from more than 25 low-mass 
NS X-ray binaries (LMXBs) can be used to put additional constraints on the high-density EoS. A pair of such QPOs is 
often seen from these systems 1 1 1]. In all currently viable models for these QPOs, the higher QPO frequency is close 
to the orbital frequency at some special radius. For such a QPO to last the required many cycles (up to ^ 1000 in some 
sources), the orbit must be outside the star. According to general relativity theory the orbit must also be outside the 
innermost stable circular orbit (ISCO). Gas orparticles inside the ISCO would spiral rapidly into the star, preventing 
the production of sharp QPOs. This implies Ill2,[l3ll that the observation of a source whose maximum QPO frequency 
is Vmax ~ Visco limits the stellar mass and radius to 

M < 2.2 MQi^ii^fl +0.75;) , < 19.5 km™i2(l +0.2j) . (1) 

- Vmax ^ ' ' Vmax \ J / v / 

The quantity j = cJ/GM^ (with J the stellar angular momentum) is the dimensionless spin parameter, which is 
typically in the range between 0.10 and 0.2 for these systems. Equation ([T]| implies that for given observed value 
Vmax the mass and radius of that source must be inside a wedge-shaped area, as shown in Fig. [T] Since the wedge 
becomes smaller for higher Vmax, the highest frequency ever observed, 1330 Hz for 4U 0614+091 [14], places the 
strongest constraint on the EoS. As can be seen from Fig.[T] the current QPO constraints do not rule out any of the EoS 
considered here. However, because higher frequencies imply smaller wedges, the future observation of a QPO with a 
frequency in the range of 1500 — 1600 Hz would rule out the stiffest of our EoS. 

If there is evidence for a particular source that a given frequency is close to the orbital frequency at the ISCO, then 
the mass is known to a good accuracy, with uncertainties arising from the spin parameter. This was first claimed for 



4U 1820-30 ifTsll . but complexities in the source phenomenology have made this controversial. More recently, careful 
analysis of Rossi X-ray Timing Explorer data for 4U 1636-536 and other sources [16] has suggested that sharp and 
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FIGURE 1. Mass-radius constraints from thermal radiation of the isolated NS RX J 1856 (grey hatched region) and from QPOs 
in the LMXBs 4U 0614+09 (wedge-like green-hatched area) and 4U 1636-536 (orange hatched region). For 4U 1636-536 a mass of 
2.0 ±0.1 Mf) is obtained so that the weak QPO constraint would exclude the NLp and NLp5 EoS, whereas the strong one renders 
only the stiffest EoS D^C, DD-F4, BBG and DBHF viable. XTE J1739-285 OH would favor a soft EoS, which is controversial, 
however (see text). 

reproducible changes in QPO properties are related to the ISCO. If so, this implies that several NS in low-mass X-ray 
binaries should have gravitational masses between 1.9M0 and possibly 2.1 M© lfT6ll . In Fig.[T]we show the estimated 
mass of2.O±O.lM0for 4U 1636-536. 

Recently, mass-radius constraints have been reported for the accreting compact stars XTE J 1739-285 1 18] and SAX 
J1808.4-3658 ifTgll (SAX J1808 for short) which are based on the identification of the burst oscillation frequency with 
the spin frequency of the compact star It is almost impossible to fulfill the constraints from RX J 1856 (and other high- 
mass candidates) and SAX J 1808 (which favors a soft EoS) simultaneously so that the status of SAX J 1808 is currently 
controversial. It is Ukely that the small radius estimate of Leahy et al. 11911 is a consequence of the underestimation of 
higher harmonics when only timing data are analyzed and not also the energy spectra. 

The nearby isolated NS RX J 1856.5-3754 (RX J 1856 for short) belongs to a group of seven objects which show 
a purely thermal spectrum in X-rays and in optical-UV. This allows the determination of R^/d, the ratio of the 
photospheric radius Roo to the distance d of the object, if the radiative properties of its photosphere are known. RX 
J 1856 is the only object of this group which has a measured distance obtained by Hubble Space Telescope (HST) 
astrometry. After the distance of 1 17 pc |20] became known several groups pointed out that the blackbody radius of 
this star is as large as 15 to 17 km. Although both the X-ray and the optical-UV spectra are extremely well represented 
by blackbody functions they require different emission areas, a smaller hot spot and a larger cooler region. The overall 
spectrum could also be fitted by blackbody emission from a surface showing a continuous temperature distribution 
between a hot pole and a cool equator, as expected for a magnetized NS. The resulting blackbody radii are 17 km 
(two blackbodies) and 16.8 km (continuous temperature distribution) ll2Tll . In this paper we adopt the result of the 
continuous temperature fit, Ro<, = 16.8 km. More recent HST observations of RX J1856 indicated larger distances of 
up to 178 pc. A distance of around 140 pc for RX J1856 is considered a conservative lower limit. For a distance of 
140 pc the corresponding radius is 17 km ll22ll . Although some questions-in particular that of the distance-are not yet 
finally settled, the recent data support an unusually large radius for RX J 1856.5-3754. 



Elliptic Flow in Heavy Ion Collisions 

The flow of matter in heavy ion collisions is directed both forward and perpendicular (transverse) to the beam axis. 
At high densities spectator nucleons may shield the transversal flow into their direction and generate an inhomogeneous 
density and thus a pressure profile in the transversal plane. This effect is commonly referred to as elliptic flow and 
depends directly on the given EoS. An analysis of these nucleon flow data, which depends essentially only on the 




FIGURE 2. Left panel: Pressure region consistent with experimental flow data in symmetric nuclear matter (dark shaded region). 
The light shaded region extrapolates this region to higher densities within an upper (UB) and lower border (LB). Right panel: 
Mass versus central density for compact star configurations, calculated using the UB and LB extrapolations of the flow constraint 
boundaries. The error bars take into account different models for the symmetry energies. The gray horizontal bars denote the 
expected mass of PSR B 15 16+02B within the la interval and a likelihood of amass below 1.72Me of 5%, resp. The vertical bars 
limit the density region covered by the flow constraint. 

isospin independent part of the EoS, was carried out in a particular model in Ref. ||23|]. In particular it was determined 
for which range of parameters of the EoS the model is still compatible with the flow data. The region thus determined is 
shown in Fig.|2]as the dark shaded region. Ref. |23] then asserts that this region limits the range of accessible pressure 
values at a given density. For our purposes we extrapolated this region by an upper (UB) and lower (LB) boundary, 
enclosing the light shade region in Fig.|2] 

Thus the area of allowed values does not represent experimental values itself, but results from transport calculations 
for the motion of nucleons in a collision ll23ll . Of course, it would be preferable to use the actual flow data as a constraint 
for each specific EoS, but a corresponding testing tool based on a transport code is not yet available. Therefore we 
adopt the results of ref. |23] as a reasonable estimate of the preferable pressure-density domain in symmetric nuclear 
matter. Its upper boundary is expected to be stable against temperature variations ll24ll . The important fact is that 
the flow constraint probes essentially only the symmetric part of the energy per nucleon Eoin). Following Ref. f23\ 
the constraint arises for a density window between 2 and 4.5 times saturation density «j. One has, however, to keep 
in mind that at high densities this constraint is based on flow data from the AGS energy regime (Ei^b ^ 0.4 — 11 
AGeV). At these energies a large amount of the initial bombarding energy is converted into new degrees of freedom, 
i.e. excitations of higher lying baryon resonances and mesons, which makes conclusions on the nuclear EoS more 
ambiguous than at low energies. Nevertheless, the analysis of |23] provides a guideline for the high density regime 
which we believe to be reasonable. The uncertainties are represented by a region of possible pressures in Fig.|2l 

An EoS which is in accordance with the flow constraint might still violate the maximum NS mass constraint. This 
is demonstrated in Fig.|2l where in particular the lower boundary of the limiting region in Fig. |2] was extrapolated 
to construct an artificial EoS (LB) in order to obtain the mass-density relation for corresponding compact star 
configurations . The resulting mass curve in Fig. [ 2] is far from reaching even typical NS mass values as the well- 
measured mass Mbi9i3+i6 = 1.4408 ±0.0003 [25]. This demonstrates that the lower bound (LB) in Fig.|2]does 
not agree with astrophysical observations, which therefore potentially narrow the band of the flow constraint. 

The maximum mass constraint demands a certain stiffness of £'o(n) in order to obtain sufficiently large maximum NS 
masses. The small influence of Es{n) on the NS mass can also be well recognized on Fig.|2] Two different symmetry 
energies, necessary to describe NS matter, were taken from the investigated EoS. They were chosen in accordance with 
the DU-constraint and gave the largest (DD-F) and smallest (D^C) contribution to the binding energy at « = 1 fm^^. 
The resulting deviations of the NS masses are shown as error bars on the curves in Fig. |2l It results in a maximum 
difference of less than 0.2 M© for the mass curves corresponding to LB. The same was done for an artificial EoS 
extrapolating the upper boundary (UB). Here, the largest error estimate of approximately 0. 1 M© is even smaller. The 



TABLE 1. A set of nuclear equations of state. The entries are: saturation density, 
Its', binding energy, ay; incompressibility, K; skewness parameter, K'; symmetry 
energy, J; symmetry energy derivative, L; Dirac effective mass, m^j. 



EoS 




ay 


K 


K' 


J 


L 


mo 




[fm-3] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[m] 


NLp 


0.1459 


-16.062 


203.3 


576.5 


30.8 


83.1 


0.603 


NLp5 


0.1459 


-16.062 


203.3 


576.5 


31.0 


92.3 


0.603 


D^C 


0.1510 


-15.981 


232.5 


-716.8 


31.9 


59.3 


0.541 


DD-F4 


0.1469 


-16.028 


220.4 


1229.2 


32.7 


58.7 


0.556 


KVR 


0.1600 


-15.800 


250.0 


528.8 


28.8 


55.8 


0.805 


KVOR 


0.1600 


-16.000 


275.0 


422.8 


32.9 


73.6 


0.800 


DBHF 


0.1810 


-16.150 


230.0 


507.9 


34.4 


69.4 


0.678 


BBG 


0.1901 


-14.692 


221.6 


-132.4 


36.3 


79.4 




DD-RH 


0.172 


-15.73 


249.0 




34.4 


90.2 


0.686 



maximum mass for UB of about 2.0 Mq again fulfills the corresponding constraint. 

NUCLEAR MATTER 

While modern nuclear matter models give a rather similar description of the saturation and sub-saturation behavior, 
they differ considerably in their extrapolations to densities above « 2ns, the regime which is relevant for NS physics 
and heavy ion collisions. This is illustrated in Table [T] where several nuclear EoS are characterized by comparing the 
parameters of an expansion around saturation as a function of the density deviation e = (n — Wj) /n, and the asymmetry 
j3 according to 

IVIost of the shown models are derived in the framework of the relativistic mean-field approach. ||26[ IztI l28l l29tl 
allowing for non-linear (NL) self-interactions of the a meson Ii30;l . For model NLp , the isovector part of the interaction 
is described entirely in terms of p meson exchange. This is different for NLp 5 where the isovector part of the 
interaction is described in terms of both p and 5-meson exchange. The latter is generally neglected in RMF models 
llsill . RMF models with density dependent input parameters (coupling constants and masses) are represented in Table 
[T]by four different models from two classes, where in the first one density dependent meson couplings are modeled 
such that several properties of finite nuclei (binding energies, charge and diffraction radii, surface thicknesses, neutron 
skin in ^"^Pb, spin-orbit splittings) can be fitted [32]. D^C additionally contains a derivative coupling which leads 
to momentum-dependent nucleon self-energies, and DD-F4 is modeled such that the flow constraint from heavy ion 
collisions is fulfilled [I23il . The second class of these models is motivated by the Brown-Rho scaling assumption [33] 
that not only the nucleon mass but also the meson masses should decrease with increasing density. In the KVR and 
KVOR models [34] these dependences are related to a nonlinear scaling function of the a-meson field such that the 
EoS of symmetric nuclear matter and pure neutron matter below four times the saturation density coincide with those 
of the Urbana-Argonne group [35]. In this way the latter approach builds a bridge between the phenomenological 
RMF models and a microscopic EoS built on realistic nucleon-nucleon forces. The RMF models are contrasted with 
several variational models for the EoS such as APR[35], WFF[ 36], FPS[37J, a relativistic Dirac-Brueckner-Hartree— 
Fock (DBHF) modellHIIij, and a nonrelativistic Brueckner-Bethe-Goldstone (BBG) model 

As shown in Fig.[3]none of these EoS predicts a maximum NS mass below the 2(7 mass limit of L68 M© for PSR 
B1516H-02B, and even at the la mass Umit of L84 M© the softest EoS NLp and NLp5 cannot be excluded. This 
underlines the important rule of the maximum NS mass as a constraint on the high density EoS. If a pulsar with a 
mass exceeding L8 — L9 M© at the 2a or even 3a level would be observed in the future, this will impose a severe 
constraint on the stiffness of the nuclear EoS. For the set of EoS tested here, only the stiffest models, i.e. D^C, DD-F4, 
BBG, and DBHF would remain viable candidates. It is noteworthy, that the flow constraint limits the stiffness of the 
EoS to be 'rather' soft, whereas the higher the constraining maximum NS mass is the 'stiffer' the EoS needs to be. 
Both together provide complementing data. 
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FIGURE 3. Mass versus central density of NS, computed for the EoS shown in Table [T] Crosses refer to the maximum-mass 
star of each sequence, filled dots mark the critical masses and central densities beyond which the direct Urea (DU) cooling process 
becomes possible. 

QUARK MATTER 
Nonperturbative QCD 

Confinement and dynamical chiral symmetry breaking (DCSB) are key emergent phenomena in QCD. Neither of 
these is apparent in QCD's Lagrangian and yet they play a dominant role in determining the observable characteristics 
of real-world QCD. A natural starting point for an exploration of DCSB is QCD's gap equation. 

S{p;^i)-^ =Z2{iy-p + iyA{pA + ipi) + m^'^) + ^{p-li), (3) 

with the renormalized self energy expressed as 

Ap-^pi) = z, j\\ii)Dpa{p-q;ii)Y7pS{q\pi)K{q.p-^P^), (4) 

where represents a translationally invariant regularization of the integral, with A the regularization mass-scale, 
Dpa{k; jJ. ) is the dressed-gluon propagator, ) is the dressed-quark-gluon vertex, and nl"'^ is the A-dependent 

current-quark bare mass. The quark-gluon-vertex and quark wave function renormalization constants, Zi 2{i^^,A^), 
depend on the renormalization point, ^, the regularization mass-scale and the gauge parameter. A nonzero chemical 
potential introduces no divergences in addition to those present in the jx — theory. Hence, the renormalization 
constants determined at /i = are completely sufficient in-medium. The gap equation's solution can assume the 
general form ll4ll] 

S{p;ix)-^ ^ iY-pA{p^,p ■ u, C^) + /74(/?4 + ili)C{p^,p ■ u, C^) +B{p^,p ■ u, C^) , (5) 

where u = (OJn). The mass function M{p^,p ■ u) — B{p^,p ■ u, i^^)/A{p^,p ■ u, i^^) is renormahzation point inde- 
pendent. There are numerous studies of this basic Dyson-Schwinger equation in-vacuum ll42[ l43ll and in-medium 
11441] ■ It is noteworthy that extant self-consistent studies of concrete models of QCD, which exhibit both confine- 
ment and DCSB in vacuum, possess coincident deconfinement and chiral symmetry restoration transitions; e.g., 
Refs. ll45l l46l WA l48i l49ll . This result appears to foUow from the crucial role played by the in-medium evolution 
of the dressed-quark self-energy in both transitions. 

Significant effort continues to be expended on determining the precise nature of the kernel of QCD's gap equation. 
A dialogue between DSE studies and results from numerical simulations of lattice-regularized QCD is providing 
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FIGURE 4. Left panel: Evolution with fl of the dimensionless quantities A(p = 0; /x), C{p = 0;fl), and M{p = 0; /x) measured in 
GeV, computed from the gap equation Right Panel: Difference in pressure between the Nambu-Goldstone and Wigner solutions 
of the gap equation. Observe, that the Wigner phase is realized for /i > 0.38 GeV, which results in a first order phase DCSB 
transition at this value. 



important information; e.g., Refs. ||53,|5l|5l|5l|5l|5l|5l|53,|5l. This body of work can be used to formulate 
reasonable Ansdtze for the dressed-gluon propagator and dressed-quark-gluon vertex in Eq. (|4]i. 

DCSB is signaled by the appearance of a gap equation solution in which the Dirac-scalar piece is nonzero in the 
chiral limit. This effect owes primarily to a dense cloud of gluons that clothes a low-momentum quark fSQ*, '6^. 
DCSB is the single most important mass generating mechanism for light-quark hadrons; e.g., one can identify it as 
being responsible for roughly 98% of a proton's mass. A system in which the gap equation's DCSB solution is the 
favored ground state is said to realize chiral symmetry in the Nambu-Goldstone mode. The antithesis is termed a 
Wigner- Weyl realization of the symmetry. Since a system's ground state is that configuration for which the pressure 
is a global maximum the pressure difference between the Nambu-Goldstone and Wigner phase reveals which of both 
is realized. A striking signature of dynamical chiral symmetry breaking is a nonzero value of the dressed-quark mass 
function M{p,pL), which entails a nonzero chiral-limit quark condensate. Naturally, a nonzero chemical potential 
destabilizes the quark condensate and, with increasing chemical potential, results in a first-order chiral symmetry 
restoring transition at jU w M(0). This is illustrated in Fig. ID Evidence suggests a coincident first-order deconfinement 
transition [5J. 



Effective Quark Matter Models 

There is no inherent obstacle that prohibits an extension of the powerful Dyson-Schwinger approach to the finite 
density domain and in fact it appears as the most promising candidate for a future systematic modeling of dense QCD 
matter in compact stars However, for explorative studies on the possibly rich phase structure of dense 

matter in QCD it is for practical reasons useful to employ effective models which preserve certain but not necessarily 
all features of QCD. Increasingly applied are models of the Nambu-Jona-Lasinio (NJL) type where chiral symmetry 
is dynamically broken in the nonperturbative vacuum and restored at finite densities and temperatures. Although the 
standard NJL model is quite simple, based on a few parameters, such as momentum cutoff and coupling constants for a 
set of interaction channels, it offers possibilities for generalizations, like a density-dependence of the cutoff parameter 
momentum-dependent formfactors |63] and an infrared cutoff to mimic confinement ll64il65ll . 

Essentially one starts from a QCD-motivated partition function 



Z(r,/i) = J ^q^qexplj'^ dx J d'^x[q{i^-mQ + ll-f)q + Xn 



(6) 
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FIGURE 5. Hybrid EoS with a low density hadronic branch described by the DBHF approach and a high density quark matter 
branch obtained from a three-flavor NJL model with color superconductivity (diquark coupling rju) and isoscalar vector meson 
coupling (t]v). 



which explicitly accounts for phases of possible interest by considering specific channels in the interaction part ^int of 
the Lagrangian. Mass gaps, pairing gaps and the EoS are then obtained from the mean-field thermodynamic potential, 
^MF — —TIuZmf- As example serves the following model ll66ll with diquark pairing in the single flavor color-spin- 
locking (CSL), two-flavor (2SC) and three-flavor color-flavor locking (CFL) channels. 




where fl — + diagy(|, — ^, — ^)/iQ + + Ag/ig is the diagonal quark chemical potential matrix and mo = 
diag f{m^,m^,m'^) the current-quark mass matrix. The relative coupling strengths of the spin-0 and spin-1 diquark 
currents, jooA = q^iCjs'^A^Aq and joi = q^iC{y\X-i + 72 A5 + YiX2)q, are essentially free parameters, but can be chosen 
according to a Fierz transformation of the one-gluon exchange interaction, Tj^o = 3/4 and Tjoi =3/8, see |67J. We 
point out, that bosonised diquark fields in QCD play an essential role in the color singlet baryon states |68]. 
The next section shows how observational constraints help to fix a model EoS obtained in this framework. 

THE QUARK HADRON PHASE TRANSITION 

A strategy to construct hybrid EoS in agreement with observational data is to adjust available coupling strengths within 
certain limits to a subset of available data. This is demonstrated in Fig|5]for a model with an interaction different 
from that in Eq.(|7]i. Here the coupling strength of the scalar diquark (770) and isoscalar vector meson (t7v) channel 
are considered as basically free parameters and have been adjusted so, that after constructing the (Maxwell-)phase 
transition from the nuclear DBHF EoS to this QM EoS the hybrid EoS fits the upper limit of the flow constraint |69]. 
Just by adjusting to this requirement the hybrid EoS fulfills all recently developed constraints from modern compact 
star observations ['2']. It is an interesting and encouraging result, that none of the present constraints on the mass-radius 
relation of NS rules out hybrid stars |70, 6]. 

As a last point we discuss ambiguities in the way the phase transition is modeled. The /3 -equilibrium condition in 
compact stars relates the chemical potentials of quarks and electrons by pLj ~ ji^ and fi^/ = flu + jJ-e- The mass difference 
between the strange and the light quark flavors > mu,m^ induces different densities of the down and strange quark. 
Charge neutrality requires a finite electron density and consequently fi^/ > jj.,,. When increasing the baryochemical 
potential, the d-quark chemical potential is the first to reach a critical value where the the partial density of free d-quarks 
becomes finite in a first order phase transition. Due to the finite value of jj.^ the u-quark chemical potential is below the 
critical value and the s-quark density is zero due to the high s-quark mass. Therefore the negative electric charge of 
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FIGURE 6. Density profiles of two stars with masses 1.4 M0 and 2.0 My,. In this model a mixed phase of d-CSL quark matter 
with nuclear matter extends up to the crust-core boundary. 



the d-quark cannot be compensated in a pure quark phase and single-flavor d-quark matter under NS conditions does 
not exit. The situation changes considering a chemical equilibrium of the type n + n p + 3d, with positive charged 
protons which compensate the charge of a single-flavor d-quark phase and results in a mixed phase of nucleons and 
down quarks once the d-quark chemical potential exceeds the critical value. This can be thought analogous to the 
dissociation of nuclear clusters in the crust of NS (neutron dripline). One interesting consequence of the effect is 
illustrated in Fig|6] Ignoring in a first approach important effects of the inhomogeneous phase mixture, e.g. Coulomb 
interactions and surface tension, the existence of such a phase suggests small fractions of quark matter even near the the 
crust-core boundary. We point out, that neither way to model the phase transition can be considered to be "correct". 
Nucleons actually consist of quarks and as expressed before it will take serious efforts to describe the QCD phase 
transition within an appropriate and fundamental approach. 
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